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Abstract 

. We continue monitoring supernova remnant (SNR) 1987A with the Chandra X- 

\ ray Observatory. As of 2004 January, bright X-ray spots in the northwest and the 

' southwest are now evident in addition to the bright eastern ring. The overah X- 

ray spectrum, Since 2002 December, can be described by a planar shock with an 
electron temperature of ~2.1 keV. The soft X-ray flux is now 8 x 10~^^ ergs cm~^ 
'■^^1 s~^, which is about five times higher than four years ago. This flux increase rate is 

■ ■ consistent with our prediction based on an exponential density distribution along 

the radius of the SNR between the HII region and the inner ring. We still have no 
^ ' direct evidence of a central point source, and place an upper limit of Lx = 1.3 x 

. 1034 

ergs s ^ on the 3—10 keV band X-ray luminosity. 



^ ' Key words: supernovae; general, supernovae; individual (SN 1987A), supernova 

■ remnants. X-rays; general 
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1 Introduction 



We present the latest results from the Chandra observations of supernova 
remnant (SNR) 1987A. As of 2004 January, we have performed a total of 
nine observations of SNR 1987A with the Advanced CCD Imaging Spectrom- 
eter (ACIS) on board Chandra (Table 1). We have presented the results from 

Email address: park@astro.psu.edu (Sangwook Park). 



Preprint submitted to Elsevier Science 



5 February 2008 



the first seven observation s in the literature ( Burrows et al et alJ 



20021 : iMichael et all liooi IPark et ID . 120041 ). Previous Chandra/ ACIS obser- 
vations of SNR 1987A have revealed a ring-like X-ray morphology with an 
asymmetric X-ray intensity distribution between the eastern and the western 
sides of the SNR. X-ray-bright spots have developed around the ring, and soft 
X-ray flux has been nonlinearly increasing. The X-ray spectrum was found 
to be thermal in origin, arising from the interac tion of the blast wave s hock 
with the circumstellar material. Previous works ( Park et al. . 20021 2004) suc- 
cessfully demonstrated that these X-ray features, when compared with the 
optical and radio emission features, were consistent with the standard phys- 
ical model: i.e., the SN blast wave shock is approaching the dense inner ring 
which was prod uced by non-spherica lly symmetric stellar winds from the mas- 
sive progenitor ( Michael et al. . 2000l ). In this model, the optically bright spots 
and the soft X-ray spots are produced by the decelerated shock entering the 
dense protrusions on the surface of the inner ring, while the radio and the 
hard X-ray emission features originate from the fast shoc k propagat i ng int o 



20021) . 



the tenuous HII region between the dense protrusions ([Park et al. 
Considering its core- collapse, Type II SN origin from a massive progenitor 
(jSonneborn. Altner. &: Kirshneii. Il987^■ as confirme d by the detection of the 



accompanying neutrino burst (jKoshiba et al.l . Il987| ) , the presence of a com- 



pact stellar remnant is predicted for SNR 1987A. There has, however, been no 
evidence for a pointlike source at the center of the remnant, which has been 
attributed to optical ly thick stellar ejecta material su rround ing t he putative 
central point source ([Burrows et al... 200Qi : IPark et al.i . 2002. |200J) . 

Table 1. Chandra Observations of SNR 1987A 



ObsID 


Date 


Instrument 


Exposure 


Source 




(Age'^) 




(ks) 


Counts 


124+1387^ 


1999-10-06 (4609) 


ACIS-S + HETG 


116 


690 


122 


2000-01-17 (4711) 


ACIS-S3 


9 


607 


1967 


2000-12-07 (5038) 


ACIS-S3 


99 


9031 


1044 


2001-04-25 (5176) 


ACIS-S3 


18 


1800 


2831 


2001-12-12 (5407) 


ACIS-S3 


49 


6226 


2832 


2002-05-15 (5561) 


ACIS-S3 


44 


6429 


3829 


2002-12-31 (5791) 


ACIS-S3 


49 


9274 


3830 


2003-7-8 (5980) 


ACIS-S3 


45 


9668 


4614 


2004-1-2 (6157) 


ACIS-S3 


47 


11856 


since the SN explosion. 



''The first observation was split into two sequences. 
As the blast wave eventually sweeps through the main body of the dense in- 
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ner ring, a dramatic increase of the X-ray flux and morphological/spectral 
changes of SNR 1987A are expected. We continue X-ray observations of SNR 
1987A with the Chandra/ ACIS in order to monitor the dynamic evolution, as 
well as to discover the central point source, of SNR 1987A. As a continuation 
of the previous works, we here briefly update the morphological and spec- 
tral evolutions of the X-ray remnant of SNR 1987A with two new Chandra 
observations. 



2 X-Ray Images 



Fig 1 displays the 0.3—8 keV band images of SNR 1987A from the Chan- 
dra/ ACIS observations. SNR 1987A has been signiflcantly brightening for the 
last ~four years. As the blast wave enters the main body of the dense inner 
ring. X-ray-bright spots flrst appeared on the eastern side of the remnant but 
have now been developing in the west. SNR 1987A has thus become closer to 
a complete ring than ever. 
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Fig. 1. Exposure-corrected, broadband (0.3—8 keV) X-ray images of SNR 1987A 
from the Chandra/ ACIS observations. Images have been deconvolved with the 
ACIS PSF and then smoothed follo wing the method described in literature 
(jBurrows et al.l . l2nnnl : IPark et a,1 .1 . [2nn2h . Darker gray-scales are higher intensities. 
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3 X-Ray Spectrum 



We perform spectral analysis of SNR 1987A following the procedure presented 
in Park et al. (2004). The X-ray spectum can be fitted with a single te mper- 
ature plane-parallel shock model ( Borkowski. Lverlv. fc Revnoldsl . 2001 ). The 
fitted foreground column and the elemental abundances are consistent within 
statistical uncertainties among the individual spectra. It is also reasonably 
expected that the foreground column and the elemental abundances would 
not significantly change over a few year-period. We thus fit the observed spec- 
tra simulaneously, assuming constant and abundances for the last four 
years, in order to provide the best constraints on the fitted column and el- 
emental abundances (Fig 2). We varied elemental abundances for N, O, Ne, 
Mg, Si, S, and Fe, while held in common for all individual spectra. We fixed 
the abundances for He and C at values appropriate for the ring abundances 
f Lundqvist &: Fransson . 1996| ). and fix ed Ca, Ar, and Ni abun dances at the 
values for the Large Magellanic Cloud ( Russell fc Dopita . 19921 ). because the 
contribution of X-ray line emission from these elemental species in the fitted 
energy band (0.4 — 5 keV) is insignificant. We varied electron temperature 
and the ionization timescale freely for individual spectra in order to trace 
any significant variations in these parameters, as perhaps expected from the 
shock evolution of this young SNR. For this purpose, we use six spectra with 
good photon statistics (>6000 counts) for statistically reliable spectral fits. 
The best-fit parameters (with 2a uncertainties) are presented in Table 2 and 
Table 3. 




channel energy (keV) 



Fig. 2. Four ACIS spectra of SNR 1987A taken between 2000 December (the lowest 
flux) and 2004 January (the highest flux), separated by ~one year from each other. 
The best-fit single-temperature plane-parallel shock model is overlaid with each 
spectrum. 
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Until 2002 December, the electron temperature of SNR 1987A h as been de 



creas ing while the volume emission measure (EM) was increasing (jPark et al. 



20041 ) ■ The latest two observations, however, suggest that the temperature 
appears to be constant since 2002 December, while the EM was continuously 
increasing (Fig 3). Implications of this recent trend of the electron temperature 
variation are not straightforward. Park et al. (2004) attributed the decrease of 
electron temperature to a physical picture where parts of the shock front enter 
the dense protrusions as the blast wave approaches the inner ring. The latest 
ACIS images, as well as the HST images, indicate that the blast wave shock 
front is now entering the dense circumstellar medium (CSM) all around the 
inner ring. The simple physical picture with small portions of the shock front 
decelerating as it approaches the dense protrusions, then, might not be ade- 
quate any longer to describe the overall X-ray spectrum. The X-ray emission 
might now be dominated by the shock entering the main body of the inner 
ring. The evolution of the overall electron temperature in this "new phase" of 
the shock is unclear. For instance, it may continue to soften, but not neces- 
sarily in the previously-reported rate, or it may establish a stable phase. The 
details of the shock evolution may be complex depending on thermal condi- 
tions of the plasma (e.g., ionization states and/or electron-ion equilibration 
behind the shock) primarily due to the radial and azimuthal density distribu- 
tion of the inner ring and its environments. We will need upcoming Chandra 
observations in order to determine whether the latest temperature variation 
means a truly stable phase or only a temporary state in the course of further, 
continuous softening of the overall X-ray spectrum. 

Table 2. Best-Fit" Elemental Abundances 



Element Abundance^ Element Abundance'' Element Abundance^ 



He 


2.57 (fixed^) 


Ne 


20+°-°2 

U-^^-0.04 


Ca 


0.34 (fixed'^') 


C 


0.09 (fixed'^) 


Mg 




Ar 


0.54 (fixed'^) 


N 


n 07+0.10 

-0.08 


Si 


n. 09+0.05 

"-'•"^^-0.05 


Fe 


n 1 5+°-°^ 
'-'•-'^^-0.01 





'-'•"-'^-0.01 


S 




Ni 


0.62 (fixed"^) 



The best-fit N^ = l.S+g;^ x lO^i cm'^. x^/i^ = 699.7/569. 
^Abundances with respect to solar. 
Lundqvist & Fransson 1996 
Russell & Dopita 1992 



Although the overall X-ray spectrum of SNR 1987A can be described by a sin- 
gle temperature plane-shock in non-equilibrium ionization (NEI), the physical 
model of the shock-inner ring interaction indicates that X-ray emission most 
likely originates from multiple velocity components of the shock in various ion- 
ization states. Considering such a multi-component shock, a two-temperature 
shock model was implemented in order to fit the observed X-ray spectrum 
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Table 3. Best-Fit Shock Parameters from A Single- Temperature ModeF'^. 



Date Electron Temperature Ionization Timescale Emission Measure 

(keV) (10^° cm-^ s) (10^^ cm'^) 



2000-12-7 




q 09+0.50 
'->-U^_0.43 


6.84l°:i 


2001-12-12 


2.48+°:i 


9 09+0.57 
^•°^-0.43 




2002-5-15 


9 9^+0.29 
^■^^-0.28 


q 9q+0.65 
"^■^"-'-0.52 


1 9 QO~l"0-'^8 


2002-12-31 


2.08l°:JI 


q 9n+0.48 
"^•^"-'-0.41 


16.68l°:?i 


2003-7-8 


9 1 1+0-14 
^•-^-^-0.13 


iq+0.42 




2004-1-2 


9 iq+0.13 

^•-"-"-•-o.is 


oq+0.45 


22.89lJ;i 



"^2/1/ = 699.7/569. 
* The best-fit N^/ = 
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Fig. 3. The electron temperature and EM variations of SNR 1987A. 

iii 



( Michael et al 



Park et al.l . 120041 ). The two-temperature shock model 



may grossly oversimplify reality, but should be a physically more realistic 
model than a single shock. In this model, the observed X-ray spectrum origi- 
nates from two characteristic components. The soft component arises from a 
slow, decelerated shock entering the dense protrusions of the inner ring. The 
hard component is produced by a fast shock propagating into the low den- 
sity HII region. The slow shock was then found to be in coUisional ionization 
equilibriu m (CIE), which w as reasonable for the shocked dense circumstellar 



material (jPark et al.Ll2004^ . Based on these previous results, we apply a two- 
component model to the latest data (Note: As discussed above, the overall 
shock phase for the latest X-ray spectrum might be different from that of 
previous observations. Nonetheless, a two-temperature model should still be 
a useful approximation for the multi-phase shock). As of 2004 January, an 
electron temperature of kT ~ 0.23 keV is implied for the slow shock in CIE 
condition, while it is kT ~ 2.2 keV for the fast shock in NEI state. The best-fit 
EMs imply electron densities Ue ~ 6300 cm~^ for the slow shock and Ue ~ 
280 cm~'^ for the fast shock. The hard to soft component EM ratio shows a 
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noticeable increase since 2002 December (Fig 4), which suggests that a signif- 
icant increase of the density in the fast shock might have begun. We note that 
the volume increase of the hard X-ray emitting plasma might also contribute 
to the observed EM ratio increase. We, however, found that the effect from 
the volume increase is negligible. This density increase is consistent with the 
blast wave, for the most parts, now entering the main body of the dense inner 
ring. 

^ 1 



o 



UJ 

CO 



LU 



0.1 



5000 6000 
Day since SN explosion 

Fig. 4. The hard to soft component EM ratios of SNR 1987A. 



4 X-Ray Lightcurve 



Fig 5 is the X-ray lightcurve of SNR 1987A. The steep, non-linear increase of 
the 0.5—2 keV band X-ray flux continues. As of 2004 January, the observed 
0.5—2 keV X-ray flux is fx ~ 8.2 x 10~^^ ergs cm~^ s~^ which is ~5 times 
higher than 1999 October. Following the same method described by Park et 
al. (2004), we fit the lightcurve assuming an exponential density distribution 
along the radius of the inner ring. The best-fit model indicates that a den- 
sity ratio between the high-density inner ring and the low-density HII region 
{interior to the inner ring) is ~18, and that a scale height of the exponential 
density distribution between inner ring and the HII region is ~0.015 pc at a 
distance of 50 kpc. The estimated density ratio is consistent with the soft to 
hard component density ratio derived from the two-component spectral fit. 



5 Point source 



After ~17 years of the SN explosion, we still find no direct evidence of a point- 
l ike source at the center of SNR 1987A . Utilizing a Monte Carlo simulation 
fjBurrows et al.Ll2000l : jPark et al.l . 120021 ) . we estimate a 90% upper limit on the 
3—8 keV band point source counts based on the latest ACIS image of SNR 
1987A. We assume a power law spectrum of a photon index F = 1.7 for the 
putative neutron star in order to estimate the X-ray flux corresponding to the 
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SN 1987A: X-ray Flux (0.5 - 2 keV) 
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Fig. 5. The X-ray lightcurve of SNR 1987A. The ROSAT data points are taken from 
Hasinger, Aschenbach, &; Triimper (1996). The dotted curve is the best-fit model 
assuming an exponential density profile. Predictions for the next two data points 
are also marked. 

derived upper limit on the point source counts. We note that the absorbing 
column toward the center of the SNR is highly uncertain depending on the 
density structure of the actual stellar ejecta prevailing there. Here we assume, 
as a first-order-estimation, a similar absorbing column Nh as measured from 
the spectral fitting of the entire SNR. We then obtain an upper limit of Lx 
= 1.3 X lO"'^ ergs s^^ on the 3—10 keV band X-ray luminosity of the cen- 
tral point source. Considering higher columns, for instance, if the absorbing 
column for the point source was an order of magnitude higher than that for 
the entire SNR, this upper limit could be a few times higher than the above 
estimate. 



References 

Borkowski, K. J., Lyerly, W. J., & Reynolds, S. P. Supernova Remnants In 
The Sedov Expansion Phase: Thermal X-Ray Emission. ApJ, 548, 820-835, 
2001 

Burrows, D. N., Michael, E., Hwang, U. et al. The X-Ray Remnant Of SN 

1987A. ApJ, 543, L149-L152, 2000 
Hasinger, G., Aschenbach, B., & Triimper, J. The X-ray lightcurve of SN 

1987A. A&A, 312, L9-L12, 1996 



8 



Koshiba, M., Thomas, M., Ryder, S., ct al. Supernova 1987A In The Large 
Magellanic Cloud. lAU Circ. 4338, 1987 

Lundqvist, P. & Fransson, C. The Line Emission From The Circumstellar Gas 
Around SN 1987A. ApJ, 464, 924-942, 1996 

Michael, E., McCray, R., Pun, C. S. J. et al. Hubble Space Telescope Spec- 
troscopy Of Spot 1 On The Circumstellar Ring Of SN 1987A. ApJ, 542, 
L53-L56, 2000 

Michael, E., Zhekov, S., McCray, R. et al. The X-Ray Spectrum Of Supernova 
Remnant 1987A. ApJ, 574, 166-178, 2002 

Park, S., Burrows, D. N., Garmire, G. P. et al. Monitoring The Evolution Of 
The X-Ray Remnant Of SN 1987A. ApJ, 567, 314-322, 2002 

Park, S., Zhekov, S. A., Burrows, D. N. et al. A Chandra View Of The Mor- 
phological And Spectral Evolution Of Supernova Remnant 1987A. ApJ, 610, 
275-284, 2004 

Russell, S. C. & Dopita, M. A. Abundances Of The Heavy Elements In The 
Magellanic Clouds. Ill - Interpretation Of Results. ApJ, 384, 508-522, 1992 

Sonneborn, G., Altner, B., & Kirshner, R. P. The Progenitor Of SN 1987A - 
Spatially Resolved Ultraviolet Spectroscopy Of The Supernova Field. ApJ, 
323, L35-L39, 1987 



9 



